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In this report, we study within the context of general relativity with one extra dimension com- 
pactified either on a circle or an orbifold, how radion fluctuations interact with metric fluctuations 
in the three non-compact directions. The background is non-singular and can either describe an 
extra dimension on its way to stabilization, or immediately before and after a series of non-singular 
bounces. We find that the metric fluctuations transfer undisturbed through the bounces or through 
the transients of the pre-stabilization epoch. Our background is obtained by considering the effects 
of a gas of massless string modes in the context of a consistent 'massless background' (or low en- 
ergy effective theory) limit of string theory. We discuss applications to various approaches to early 
universe cosmology, including the ekpyrotic/cyclic universe scenario and string gas cosmology. 
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I. INTRODUCTION 

The idea that the universe is comprised of any num- 
ber of extra dimensions, in addition to the three non- 
compact spatial dimensions that we observe, goes back 
all the way to the work of Kaluza and Klein in the 1920's 
P, @ , where extra dimensions were first proposed as a 
manner in which one could unify particle interactions 
with gravity. Since then, theories with extra dimensions 
have evolved from being a mathematical curiosity, albeit 
one with remarkable consequences for particle physics , 
to a theoretical necessity with the advent of string the- 
ory. Since the famous anomalies cancellation calculation 
of Green and Schwarz the challenge has remained 
for cosmologists and phenomenologists alike to account 
for the presence of extra dimensions in physics. From 
a pessimist's point of view, one might view these extra 
dimensions as an unwanted surplus that has to be swept 
under the rug at presently accessible energies. In its most 
unimaginative form, one could view this as the gist of the 
moduli problem. From an optimist's point of view how- 
ever, one could view the requirement of extra dimensions 
for the consistency of string theory as a new window of 
opportunity, through which one might be able to resolve 
some of the outstanding problems of particle physics and 
cosmology. 

Taken as an essential part of the radically new pic- 
ture of space-time that string theory offers us, extra di- 
mensions, with strings and branes of various dimensions 
propagating through them, could offer us a potentially 
non-anthropic explanation for the observed dimensional- 
ity of space-time In addition to offering us poten- 
tially non-singular cosmologies [a, Ll 13 > they could very 
well be at the heart of the resolution of the dark mat- 
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ter problem [1 H, OH El . Recently, several alternative 
stringy cosmologies have emerged as possible successors 
to the standard big bang model of the universe, such 
as the pre-bi g bang, model 0, brane/string gas cosmol- 
ogy J5 |. tq . [ Ta. , a nd the cyclic/ekpyrotic scenarios 
UJA llo. fl7l Ha . fl9t l20j | , all of which have the ultimate 
aim of becoming complete and testable models of the 
early universe. 

A central question that arises in determining whether 
or not these models of the early universe reproduce ob- 
servations is, how metric fluctuations in the presence of 
extra dimensions evolve in the backgrounds that these 
cosmologies propose. The importance of understanding 
this question cannot be understated in the context of the 
abundance of available experimental data, against which 
we must compare our eventual predictions. For instance, 
in the ekpyrotic/cyclic scenario the effects of branes col- 
liding in a 5-dimensional bulk are explored [60| 

Key to the success of this program is being able to 
follow the evolution of metric fluctuations through the 
bounces of the extra dimension. The eventual goal of 
this is to be able to explain the observed inhomogencities 
of the cosmic microwave background as being seeded by 
radion fluctuations generated through the motion of the 
branes towards each other, rather than through the quan- 
tum fluctuations of some putative inflaton field. 

Similarly, in the context of the pre-big bang scenario, 
where an expanding and a contracting phase of the uni- 
verse's evolution are naturally related to each other by 
the scale factor duality symmetry inherent to low energy 
effective string theory, a universe undergoing a big bang- 
big crunch transition is naturally implemented. In this 
context, the evolution of metric fluctuations through the 
bounce becomes a question of central importance in un- 
derstanding the way the big bang phase of our universe's 
evolution happened the way it did. As it turns out, al- 
though at first seemingly unrelated, an associated and as 
yet unexamined issue is how any tentative solution of the 
moduli problem will modify the spectrum of the metric 
fluctuations which we observe. Given that the universe is 
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likely to contain many extra dimensions whose shape and 
volume moduli are dynamical in the early universe, one 
might expect that the dynamics of these moduli fields 
will lead to a highly nontrivial modification of any initial 
pre-stabilization spectrum. 

It is the goal of this report to show that, in the context 
of a 5-d universe where an extra dimension undergoes 
non-singular bounces, or is en route to stabilization, the 
final spectrum for the Bardeen potential corresponding 
to (long wavelength) fluctuations of the scale factor for 
the non-compact dimensions, is identical to the initial 
spectrum of the five-dimensional fluctuations (in a sense 
which will be specified later). Specifically, if we had an 
initial scale invariant spectrum for the five-dimensional 
fluctuations, we would end up with a scale invariant spec- 
trum for the Bardeen potential after the bounce, or af- 
ter the transients of the stabilization had settled down. 
Furthermore, we show that the variable corresponding 
to metric fluctuations of the compact dimensions decays, 
which corresponds to the stability of the radion degree of 
freedom to fluctuations in our non-singular setup. 

One application of our work concerns the evolution of 
fluctuations in ekpyrotic/cyclic type models. In the con- 
text of a four space-time dimensional effective field the- 
ory toy model of this scenario, the dynamics of perturba- 
tions has been investigated in detail. The initial analyses 
[l5ll22| yielded the result that a scale- invariant spectrum 
before the bounce transfers to a scale-invariant spectrum 
after the bounce. These analyses, however, were ques- 
tioned inHHHIllllllllllil. A serious complicat- 
ing factor turned out to be the fact that the proposed 
background evolution was singular, thus requiring the 
use of "matching conditions" (such as those derived in 
|29j in the context of an expanding cosmology undergoing 
a sudden phase transition) to compute the post-bounce 
spectrum of fluctuations. As discussed in [3(j], there is a 
very sensitive dependence of the final result on the choice 
of the matching surface. Independently, there has been 
recent work on the evolution of fluctuations through a 
nonsingular bounce in four space-time dimensional cos- 
mologies, in which the bounce is constructed by adding 
extra terms to the standard Lagrangian. These analyses 
HH H3, H3, [S3, HE HE E3, Ea yield results show- 
ing a sensitive dependence on the nature of the bounce. 
The bottom line of this work is that a correct analysis of 
fluctuations in the ekpyrotic/cyclic scenario needs to be 
done in a five space-time dimensional context, a context 
in which the nature of the bounce is unambiguous. A first 
important step in this direction was taken in |4d|. con- 
firming the result that a scale-invariant spectrum passes 
through the bounce (in the five-dimensional context, the 
bounce means that the radius of the extra spatial dimen- 
sion bounces, not that the four space-time dimensional 
scale factor bounces) without change in the spectral in- 
dex. However, the analysis of 01 was done in the context 
of a singular background and assuming specific match- 
ing conditions for fluctuations applied at a point when 
the perturbations in fact blow up. Thus, the results are 



open to doubt. In this work, we study the transfer of 
fluctuations through a cosmology in which two boundary 
branes approach each other and bounce without encoun- 
tering a singularity (see also [2(]] for previous work done 
in the context of a particular nonsingular ekpyrotic-type 
bounce proposed in 41]). 

We begin by introducing a non-singular bouncing 
model of a 5-dimensional universe where one dimension 
is compactified on a circle or an orbifold. As we shall 
see shortly in detail, our non-singular background is ob- 
tained by considering the effects of gas composed of mass- 
less string modes on the dynamics of space-time, in the 
context of a consistent low energy effective theory limit, 
or 'massless background', of string theory. We then set 
up the framework for studying cosmological perturba- 
tions in this model and study how these transfer through 
the various bounces that the extra dimension undertakes. 
Although our setup is seemingly specific to string gas cos- 
mology, the essence of our framework is that we have a 
non-singular bounce/stabilization mechanism that is af- 
fected by degrees of freedom that become massless at a 
certain point. Hence, the hope is that the results ob- 
tained here can be generalized to other settings, an issue 
we will discuss in detail when we consider applications 
to different approaches to stringy cosmology. We now 
commence our paper with a few preliminaries. 



II. THE BACKGROUND 
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Consider a five dimensional space-time with the topol 
ogy of R 4 x S 1 , described by the metric [6l| 

diag(-l, a 2 (t), a 2 (i), a 2 (t),b 2 (t)) , 
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from which we derive the following components of the 
Einstein tensor: 
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Here, the indices i and j run over the three large spatial 
dimensions and y denotes the extra dimension. The form 
of the energy-momentum tensor which will couple to the 
Einstein tensor is given by 
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We can recast the Einstein equations G A B 
where G is the five-dimensional Newton's constant, in the 
form: 
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where H = a/a and TL = b/b. The energy-momentum 
tensor for a string gas in this toroidally compactified 
background was derived in [42j | to which we refer the 
reader if any of what follows is unfamiliar. It was found 
to be 
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where g denotes the determinant of the metric, b := 
bj\fa! (with 2na' being the inverse string tension) and 
the subscripts n, w refer to the momentum and wind- 
ing quantum numbers of a closed string along the 5th 
dimension respectively, on which all of the subscripted 
quantities will depend. The density fiQ is the number 
density of the string gas with the metric dependence fac- 
tored out (p(t) = po/y/—g), while e n , w is the energy of a 
single closed string in this background 
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where p n . c . denotes the center of mass momentum along 
the three non-compact directions, and N is the number 
of right-movingstring vibrational modes. For reasons 
made clear in [l^] and 0], string states that are mass- 
less at special symmetry points (i.e. the self-dual radius) 
should play a very distinguished role in any string gas 
cosmology (see also EH)- There are several rea- 

sons for this, among which are prominently the desire to 
obtain a viable late time phenomenology and a robust 
stabilization mechanism for the radion. However, the 
main motivation is that the low energy effective theory 
limits of string theory (such as general relativity and dila- 
ton gravity) are also consistent backgrounds on which we 
can study the propagation of massless (and only mass- 
less) strings. For these reasons we do not further justify 
the focus on massless string modes, and refer the reader 
to 0> HE E3 f° r more details. For a general review of 
stri ng gas cosmology, we refer the interested reader to 

We wish to emphasize that the framework within which 
we chose to work assumes nothing other than the fact 
that nature is described by a string theory at high en- 
ergies, and that the degrees of freedom that are likeliest 
to be excited (namely massless string modes) will be ex- 
cited. The two particular geometries that we consider for 
the extra dimensions (i.e. a toroidal geometry in this sec- 
tion and an orbifold geometry in the next section), are 
necessitated by the fact that they are the unique com- 
pact backgrounds on which one stays within the approx- 
imation of the low energy effective limit of string the- 
ory, when considering extra dimensions that are similar 
in size to the string scale. This is because in general, 



were the metric to depend on the extra dimension (as is 
the case for warped extra dimensions), one will introduce 
curvatures that approach the string scale when the ex- 
tra dimensions themselves approach the string scale, and 
hence the low energy approximation will break down. We 
discuss this point further in the next section. 

The particular states we are interested in are those 
for which n — —w = ±1,N = 1. That these states 
are massless at the self-dual radius (6 = *Ja!) is easily 
checked from (|12|l . These states imply that {7J becomes 
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We will see that this background stabilizes around the 
self-dual radius (b = 1), whilst the non-compact direc- 
tions persistently expand as they would in a radiation 
dominated universe. Before we get to this, however, we 
wish to discuss a generalization of this background to 
the situation where the extra dimension is compactified 
on an orbifold. 



Extension to 



The consistency of the framework that we propose, 
namely that we remain within the low energy effective 
theory (or massless background) limit of string theory, 
requires the background we consider to satisfy the con- 
dition 
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That is, the Ricci scalar should be considerably bounded 
from above by the string tension. This requirement 
translates into the statement that the metric should not 
change very much on the string scale, which ensures that 
this background remains a consistent background for the 
propagation of massless string modes. Conversely, were 
one to consider backgrounds which do not satisfy i|14[l . 
not only would the approximations made in deriving 
the low energy limit of string theory break down, one 
would also expect massive modes to be created H3 . For 
our purposes, the requirement of l|14|) means that were 
we to consider compactifications of the extra dimension 
on scales comparable to the self-dual radius, then there 
should be absolutely no dependence of the metric along 
the compactified direction. This is a remarkable feature 
of using a toroidal compactification: one can in fact study 
string scale processes without invalidating the approxi- 
mations inherent in the low energy limit of string theory. 
It was shown in |43j and |47| that the background we are 
about to derive, does in fact satisfy (|14f) throughout its 
dynamics. 

Returning to the problem at hand, we see that since 
we are not allowed to consider metric dependencies along 
the extra dimension (by homogeneity and isotropy of the 
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non-compact dimensions, the metric can then only de- 
pend on time), the components of the Einstein tensor 
®-® are unchanged after orbifolding the extra dimen- 
sion. The only place where orbifolding might make a dif- 
ference is in the energy momentum tensor of the string 
gas. However, because the massless states that we have 
focused on are in the so-called untwisted sector [48| (in 
general, twisted states are localized at the orbifold fixed 
points), it turns out that there is no difference in the 
energy-momentum tensor either <j§)l- (|ll[) . This is a con- 
sequence of the 'inheritance' principle of orbifold theories 
|48|. and permits us to use the framework just presented 
in either situation, provided we only use the modes that 
we have indicated and we do not consider any variations 
in any metric quantities along the extra dimension. How- 
ever, there is a caveat to this in that two orbifold fixed 
planes are present, which might have matter confined to 
live on them. The requirement not to induce any met- 
ric variations along the extra dimension translates into, 
via the Israel junction conditions |49l |5C[. the condition 
(among others) 



da 
(Iff 



(15) 



where [da/dy] is the jump of the derivative of the scale 
factor a along the extra dimension, evaluated around ei- 
ther of the orbifold fixed planes. Hence, we see that if 
we require there to be no dependence of the metric along 
the extra dimensions, then any matter localized on the 
branes must be sufficiently dilute to render the branes 
to behave as test branes. Bearing all of this in mind, 
we can now proceed to derive the background solution 
within our framework. 



III. BACKGROUND SOLUTION IN 5D 

A. The three large dimensions 

We first note from Q that any form of matter which 
satisfies a radiative equation of state (p = p/3, r = 0) for 
all times (as opposed to our wound string states which 
become massless, hence radiative, only at the self-dual 
radius) drops out of the driving term in the equation 
of motion for b (01. Hence, we can safely consider a 
situation where we have a radiation gas in addition to 
our gas of wound strings, which dominates the evolution 
of the Hubble factor H for the non-compact dimensions 
through (jSJ) after the extra dimension has stabilized, or 
is close to stabilization. In fact, massless unwound closed 
string states (gravitons: n = w = 0, N = 1) can provide 
exactly such a bath. Henceforth, it is natural to include 
such a bath in our setup. We do not consider this matter 
further, simply taking as a given in what follows that the 
scale factor a for the non-compact dimensions expands at 
the background level as it would in a radiation dominated 
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FIG. 1: The analytic approximation 1231 (solid line) is com- 
pared with the numeric solution of l|2()jl (circles). 



universe, that is 
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B. The extra dimension 



Let us go back to the equation of motion l|13|) for 6, 
which reads 
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where p n . c . =■ q/a with q = const is the center of mass 
momentum along the three large dimensions. Introduc- 
ing the characteristic time scale 
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and defining the dimensionless time variable 




such that a = agr), we can simplify (|17fl to 



b" + -V 

v 



-J- + h - + 2 

i>2 ' 3 ' 3 



(18) 



(19) 



(20) 



where a prime denotes a derivative with respect to ry and 
we also introduced the free dimensionless parameter 



I 



fa' q 



(21) 



This equation can easily be integrated numerically, but 
it will turn out to be useful to have a simple analytic 
approximation. 

If we are close to the self-dual radius, that is b = 1 + e 
with e -C 1, we can expand the driving term in (|20l) so 
that 
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results. The general solution is given by 
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e(rj) = — (vising) +Bcos(rj)) , 
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with A and B constants. This solution may either de- 
scribe the stabilization of the extra dimension or a series 
of non-singular bounces j^]. See Fig for a compari- 
son of the numeric solution of l|20() with the analytic one. 
Since they are barely discernable we will use the analytic 
approximation in the following and switch freely between 
bounce / stabilization-language. 



IV. SCALAR PERTURBATIONS IN 5D 

We will first perturb the metric, focusing only on scalar 
metric perturbations. The most general scalar perturba- 
tion in "generalized" longitudinal gauge [2(J, I&J |52| can 
be characterized by four scalar functions, $, ^, W and 
r . These functions can be viewed as a basis of gauge 
invariant variables. In this gauge, the metric is given by 

ds 2 = -(1 - 2$(i,x))(ft 2 +a(£) 2 (l - 2*(<,x))dx 2 

+6(t) 2 (l + 2T(t, x)) dy 2 - 2W(t, x) dt dy , (24) 



where the signs in front of the perturbations are a mere 
convention |63|. For reasons discussed in sections UTI and 
III Al we do not consider any y dependence, that is we 
assume homogeneity in the extra dimension (see also [l4T | 
for a discussion of this and other assumptions of SGC). 
The corresponding Einstein tensor can be computed to 
be [23, E3 (in the following A denotes the Laplacian of 
the three large spatial dimensions): 
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One can check that the equations of motion involving 
W decouple from the other ones for the matter content 
iJHJl we consider. Since W would appear only squared in 
a four dimensional effective theory, we will not need to 
compute it at all j64[. 

To write down the perturbed Einstein equations |65j . 
we also need the perturbed energy-momentum tensor 
ST -g . It will include the thermal bath of radiation (de- 
noted by the subscript r), and the stringy matter sources 
denoted by a tilde. To be specific, we have 

/ 5p( r ) -(P(r) +P(r))V t \ 

(S (r) T A B ) = ( P{r} + P(r } )V t S P{r) S) , 
\ / 



(where pt r \ and pt r \ are the radiation energy density and 
pressure, respectively, and V is the radiation three veloc- 
ity potential) and 

/ Sp ~(P(r) +P(r))V,t \ 

(Sf A B ) = (p (p) + P{r ))% SpS) , 

V Sf J 

where p and p are the string gas energy density and pres- 
sure, respectively, and V is the string gas three veloc- 
ity potential. Note that anisotropic stress does not fea- 
ture in our setup (see appendix 0, but we will keep V 
around for the time being. We will also focus on adi- 
abatic perturbations of the radiation fluid only, that is 
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5pr r ) = (5/9( r )/3, even though two ideal fluids are present 
so that iso-curvature perturbations could arise. The rea- 
son for neglecting those is simplicity. The xi — Xj (i ^ j) 
Einstein equations yield immediately 



r = * 

and, after introducing 

£ := 
N := 



# - $ 
Mo 



(32) 

(33) 
(34) 



we infer from |([J}-||TTJ, that the perturbed energy mo- 
mentum tensor is of the form 
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where the background quantities are given by 
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In performing this calculation, we started with the source 
action for this string gas |42j, and obtained the per- 
turbed energy momentum tensor through direct calcu- 
lation. The string gas energy momentum tensor only 
depends on the metric, the number density of strings po, 
and the center of mass momentum of these strings. How- 
ever, the latter quantity is not perturbed, because for 
long wavelength perturbations (compared to the string 
scale), which we are restricted to if we are to remain 
within the limits of low energy effective theory, the center 
of mass motion of a string propagating on a perturbed 
spacetime is unaffected to first order. Hence the only 
quantities left to perturb are the metric dependencies, 
and the number density of the string gas. 

Next, we can write the perturbed zero component of 
the conservation equation V a Tab = for the stringy 
matter as 



= 5p + 3H(Sp + Sp) + H(6p + 8r) 



(41) 



- Utf + P) - ~f(p - 3p + 2f ) + -^(p + p) . 



For simplicity we will set V — from now on, that is 
we neglect the scalar velocity potential of the string gas, 



consistent with the arguments of appendix^] Combining 
the diagonal Einstein equations yields 
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The last three equations (|41|l - l|43|l are the dynamical 
ones for the perturbation variables T, £ and N . The it—t) 
equation gives the radiation fluid perturbation Sp( r ) and 
the (xi — t) equation gives the scalar velocity potential V 
of the radiation fluid in terms of the other variables. 

These equations simplify if we write them in terms of 77 
defined in 119|) and make use of the background equation 
a = aQf]. Introducing h :— b' jb we get 
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Here, we made the transition to Fourier space and we 
performed the rescaling a* := a/ao, k* := kto/ao. All 
other starred quantities are defined via 
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The two dimensionless constants left are 
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A. Analytic late time solution 

After h = ti = got approached (that is after the extra 
dimension got stabilized or, in the language of the ekpy- 
rotic/cyclic scenario, the two approaching branes came 
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to a halt) , we can simplify the equations of the previous 
section to 
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where we introduced the constant 
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We note that the source term for vanishes because of 
the background equation of state for the gas of massless 
string modes. Also, the equations for and Tk decouple 
at late times, just as they should do. The solutions are 
given by 

T fc = - (Ci cos(w r ?7) + C 2 smfuirr])) , (54) 
60 = ^3 (c 3 (k*r]cos(uj ( ri) - 2\/3sin(^?7)) (55) 
+C4(fc*?ysin(^77) + 2Vicos(uj^r]))] , 
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when a long wavelength mode enters the Hubble radius 
again at k*rj r = 2. If we Taylor expand Ij55(l for small 
k* we see that there is an approximately constant mode 
present for On the other hand, from (|54|) is os- 
cillating and decaying cx I/77. Therefore, whatever the 
spectrum for was at the initial time, it should per- 
sist till re-entry and determine the spectrum for $k be- 
cause, by neglecting r&, we have $t w — ~ — 6»/2- 
This is the same result one would conclude in a simple 
four dimensional universe, dominated by radiation. To 
be specific, we can approximate 



\MVr)\ 



(kjrji) j_ r fc fa) cos(fc?7 r /2 + 0) 
2 



2r]r/r]i cos{kiji/2 + P) 



,(59) 
(60) 



where (3 is some irrelevant phase. Of course this holds 
true only, if the transient stabilizing epoch leaves no 
strong imprints onto the spectrum; hence, we will ex- 
amine this crucial issue in the next section. 

Before we do that, let us have a quick look at the 
physics after a mode entered the Hubble radius: will 



start to oscillate with a decaying amplitude cx I/77 2 . If 
one waits long enough, Tk will catch up and become 
larger than Now w VPfe « Tk/2 and thus Tk will 
determine the spectrum, because for r\ r; r 



I** fa) I 



rfcfa) cos(fc7?/2 + (3) 



Zv/Vi cos(%/2 + /3) 



(61) 



There is of course a transition region where neither ap- 
proximation holds. Nevertheless, this only occurs after 
Hubble radius crossing and shall not concern us in the 
following. 



B. Numerical solution and contact with late times 

The analytic solution of the previous section is valid 
after b got stabilized (or the branes came to a halt). How- 
ever, the main question we would like to address is the 
following: what effect on the perturbations, if any, does 
the process of trapping the modulus b have? In order to 
address this question we have to solve (|44H - H46|) numeri- 
cally. As we shall see, the late time solution is approached 
rapidly and long wavelength perturbations are not af- 
fected in any significant way. Hence, the naive intuition 
that the value of $^ for modes outside the Hubble radius 
stay frozen is indeed valid, with important consequences 
for cosmological model building, be it in the framework 
of the ekpyrotic/cyclic scenario, or within string gas cos- 
mology (SGC) - we will focus on concrete models in the 
next section. 

Firstly, let us specify constants and initial conditions: 
the only free parameter in our setup is q* defined in l|49|l . 
describing the initial momentum of the string gas in the 
three large dimensions. Being a dimensionless parameter, 
we choose the most natural value 



q* := 1 . 



(62) 



Our results are not very sensitive to its exact value [66j. 
We start close to the self-dual radius with some arbitrary 
velocity (once again, its exact value does not affect our 
results - see FigEl, that is 



Hih) = 1.1, 
V{ m ) = -0.1 



(63) 
(64) 



Next, we have to specify initial conditions for ^fc,Ffc 
and Nk at 77$ = 1 |67| . Since we did not compute how the 
universe came close to the self-dual radius (e.g. within 
the ekpyrotic/cyclic setup or a SGC setup), we have no 
way of deriving those. However, we can make an ed- 
ucated guess: if the universe underwent some phase of 
inflation in the three dimension before (we will make this 
point more concrete in the next section), one should ex- 
pect similar values for all metric perturbations; hence we 
will set 



£k(r]i) = Tk{r)i) ='■ e k 



(65) 
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FIG. 2: log ( | <&fe ) | is plotted for different values of k* , with 
the initial conditions given in section II V Bl Black: analytic 
solution of l|59|l : Grey (bending curve): numerical solution of 



with some small e^. For simplicity, we will use e k = 0.01 
for all k. Furthermore, since the long wavelength modes 
we are interested in should be frozen once they cross the 
Hubble radius, we set 



(66) 



The initial string density perturbation Nk(j]i) will be set 
to zero, because these massless modes just got produced. 
Naturally, N k will get sourced by Tfc and £fc. No con- 
clusion in the following is sensitive to the chosen initial 
conditions, hence we chose the most simple ones above 
(other ones were of course also tested). 

The only thing left to specify is the wave-number k*: 
We are interested in long wavelength modes, that is 
modes outside the Hubble radius with small k*. In ad- 
dition, we shall require k* < k, so that we can compare 
our results with l|59|) . With all initial conditions specified, 
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FIG. 3: The spectrum of (a) £*, (b) Y k and (c) $ fc = 
(Tk — £*0/2 is evaluated at Hubble radius crossing r\ r = 2/k* , 
with the initial conditions given in section HVBl If the long 
wavelength regime k* <C 1 is approached, all oscillations are 
damped away such that only the constant mode of £k survives, 
as expected. 



we can use the analytic solution ill'MI) to solve 14411 - 146|) 
numerically. 

Let us first compare &k(v) = {^k(v) ~ Cfc C 7 ?)) /2 for 
different values of fc*, Fig. |21 (a)-(c): There is plotted 
both numerically (grey) and analytic (black). We see 
that the analytic late time solution is approached pretty 
fast - in fact, the small difference at the beginning is not 
visible in this plot, but only if is plotted alone as in 
Fig. [S] (a) . The visible decaying oscillation of is the 
decaying mode of form (|54|l with frequency lut, which 
is plotted in Fig. [S] (b) (compare with Fig. [21 (b)). In 
addition, there are strongly damped oscillations on top 
of the constant mode of that carry the same frequency 
as b(i]), Fig. 21 As mentioned before, these are not visible 
in Fig.[2but in Fig. 0(a) where is plotted alone. These 
oscillations are the impact of trapping the modulus 6, an 
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FIG. 4: (a) and (b): b over fc*»7 is plot ted fo r k* = 0.005, with 
the initial conditions given in section H V Bl (c): is plotted 
over the same time range as b; Note how the oscillations in 
the background scale factor source transient oscillations in 



impact that can safely be ignored in the long wavelength 
regime: if one decreases k*, all of these effects become 
less pronounced. 

An other feature visible in Fig. [3 is a bending of the 
lower curve in comparison to the analytic solution (|59|l . 
This bending is the beginning of an oscillation of and 
hence an expected feature, since the constant mode is 
actually the beginning of an oscillation with a very small 
frequency due to u>£ . All features are not sensitive to the 
initial values we choose for Tk and 

Last but not least, Nk is plotted in Fig. 0(c): it gets 
sourced quickly by Tk and and oscillates with a nearly 
constant amplitude thereafter, with all frequencies enter- 
ing. However, there is no mentionable back-reaction of 
Nk on the metric perturbations. Hence, one could give 
Nk a non zero initial value without changing the evolu- 
tion of Tk and in any significant way. 



FIG. 5: The perturbation variables Tk and Nk plotted for 
k* — 0.005, with the initial conditions given in section II V Bl 



V. CONSEQUENCES FOR 
EKPYROTIC/CYCLIC MODELS AND SGC 

In the end, we are interested in the spectrum of $k- 
So far, there seems to be little to no impact of the tran- 
sient stabilizing era of b onto Tk and This translates 
directly to the spectrum, where should approach a 
horizontal line for small k* . This is indeed the case, as 
can be seen in Fig. (c), where the spectrum is evalu- 
ated at horizon crossing: the oscillations that are present 
for relatively large k* get damped once the long wave- 
length regime is approached. This is expected, because 
long wavelength modes enter the Hubble radius later and 
henceforth, the decaying modes responsible for the oscil- 
lations in the spectrum get damped more. The frequency 
in the spectrum of <&k is proportional to 1/fc*, in accor- 
dance with r\ r ~ 1/fc* and equation (|59p. 

Note that no shift of the spectral index is induced, 
leaving the overall index unaltered. The oscillations in 
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FIG. 6: log(|$|)| is plotted for k* = 0.005. This is the same 
plot as in Fig|3](b), but with different initial conditions for b: 
b(r/i) — 1 and b 1 (rji) — 0.1 are used instead of IfiM^ and l|64l . 
Note that there is no discernable difference between this plot 
and FigH(b). 



£fc in Fig.|5Ja) are the main impact of the transient era of 
a dynamical b. They are clearly negligible, because the 
tiny amplitude gets even smaller with decreasing k*. 

To summarize, the main conclusion of our numeri- 
cal study is that long wavelength metric perturbations 
quickly approach their asymptotic solution, leaving no 
trace of the many bounces b experiences on its way to 
stabilization. More specifically, we have shown that the 
initial spectrum of the non-decaying mode of is pre- 
served, i.e. equals the spectrum of at late times. At 
late times, the spectrum of equals the spectrum of the 
Bardeen potential the potential which is relevant for 
late time observations. Speaking more loosely, we have 
shown that the initial pre-bounce spectrum of bulk per- 
turbations is transferred to the late time cosmological 
perturbations without a change in the spectral index. 

We saw in the previous section that there is no signifi- 
cant effect of trapping the modulus b on the spectrum of 
long wavelength perturbations in the sense that the final 
spectrum of the Bardeen potential $ equals the initial 
value of the non-decaying mode of the bulk perturbation 
£. This comes about since the value of £ for small k* 
modes which have a wavelength larger than the Hubble 
radius remains frozen during the phase of the trapping 
of 6, and since the final spectrum of £ coincides with the 
spectrum of 

Even though our result agrees with a "naive" intu- 
ition coming from the analysis of fluctuations in expand- 
ing four space-time dimensions, it could be viewed as 
an unexpected result: The "naive" intuition turned out 
to be wrong in the case of bouncing cosmologies in four 
space-time dimensions, for example in pre-big-bang mod- 
els where it was shown 55] that the growing mode of $ 
in the contracting phase couples almost exclusively to 
the decaying mode of $ in the expanding phase, leading 
to the result that there is a large change in the spec- 
tral index of the dominant mode between the contract- 
ing and the expanding phase. Other analyses performed 
within the context of four-dimensional general relativity 



(see [36|, |22l |3g, |39| and references therein) yielded a sim- 
ilar result. It is now also generally accepted that the four 
space-time dimensional toy models which were proposed 
to describe the ekpyrotic/cyclic models have the same 
feature, namely that the dominant mode of $ in the con- 
tracting phase matches predominantly to the decaying 
mode of $ in the expanding phase (see e.g. 1241 1. Our 
analysis, however, confirms the analysis of |40j ]. which 
showed, in the context of a singular bounce, that in the 
case of a bounce of boundary branes in five space-time 
dimensions the spectrum of <E> could be preserved. 

Our result is an important step in constructing a vi- 
able alternative to standard scalar field-driven inflation- 
ary models (next to the remarkable but highly nontriv- 
ial KKLMMT construction, [5|J and follow-up papers) 
in the context of cosmologies with extra dimensions. 
We will now outline two possible proposals which seem 
promising to us in that regard: 

Firstly, one could envisage a modification of the ekpy- 
rotic/cyclic scenario. In the ekpyrotic proposal 0, a 
test brane moves slowly through the bulk towards the 
boundary brane of an orbifold (the Horava-Witten 
setup). In the cyclic scenario [l6|]> it is the two bound- 
ary branes that approach each other. In either case, a 
scale-invariant spectrum of perturbations can be gener- 
ated during the contracting phase given a suitable po- 
tential for the modulus field which describes the inter- 
brane distance j^. During the collision of the branes 
(a singular event in the ekpyrotic/cyclic scenarios) a hot 
universe is supposed to emerge on the boundary brane 
we live on today. Our background construction can be 
viewed as a regularized version of the ekpyrotic/cyclic 
scenarios (ekpyrotic in the context of a single bounce, 
cyclic in the case we consider the evolution until b has 
stabilized). Our work shows that in the context of such a 
regularized scenario, an initial scale-invariant spectrum 
could pass through the bounces and thus survive from 
the initial contracting phase to the final phase of expan- 
sion of our three spatial dimensions (this scenario will be 
developed further in a follow-up paper). In our scenario, 
the branes do not actually hit each other, but come to a 
halt as a consequence of the appearance of new massless 
modes that get produced explosively at a certain brane- 
separation (that is at a certain value of b). Our results 
show that such dynamics does not spoil a scale invariant 
spectrum generated during the initial brane movement. 
In this framework, one can still reheat the Universe, since 
the "kinetic energy" of the branes will get transferred 
to the stabilizing massless modes and other light modes 
which act like radiation on the brane (corresponding to 
the radiation bath that we already included in our setup). 

Secondly, an incorporation of inflation into 
brane/string gas cosmology will rely on the stabi- 
lization mechanism described in this article: if any 
mechanism of inflating three dimensions is found, the 
internal dimensions will most likely have to deflate in 
some way. For example, one could employ the idea of 
anisotropic inflation (studied in the context of vacuum 
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solutions of higher dimensional general relativity in |58( 
and Once the internal dimensions get close to the 

self-dual radius, certain string modes become massless 
and get produced explosively along with a radiation 
bath, in close analogy to reheating after standard scalar 
field-driven inflation. These modes will then stabilize the 
internal dimensions, curing the graceful exit problem of 
anisotropic inflation, while the spectrum of fluctuations 
produced during inflation in the metric degrees of 
freedom remains unaltered. A concrete realization of 
this proposal is in preparation by the authors of this 
article. 



VI. CONCLUSIONS 

In this report, we studied how radion, matter and met- 
ric fluctuations interact in a universe that exhibits a tran- 
sient stabilizing epoch of its extra dimension, which could 
also be viewed as a series of bounces of the extra dimen- 
sion. We were primarily interested in the imprints of this 
epoch on the spectrum of the Bardeen potentials, moti- 
vated by the hope of combining an earlier phase, gener- 
ating a scale invariant spectrum (e.g. in the framework 
of a modified ekpyrotic scenario, or via an incorporation 
of inflation within string gas cosmology), with a success- 
ful late time stabilizing mechanism, provided by a gas of 
massless string states. 

We found that the spectrum of long wavelength per- 
turbations remains unaltered by this epoch, which has 
important consequences for various approaches to stringy 
models of the early universe. In particular, specific re- 
alizations of the ekpyrotic/cyclic scenario or an incor- 
poration of inflation into string gas cosmology become 
possible - we provided two proposals in this report, but 
many variations are viable. 
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APPENDIX A: ANISOTROPIC STRESS 

In this appendix, we will provide arguments as to why 
no anisotropic stress arises in the setup of > which we 
examined at the perturbative level in this article. In the 
following a,b, . . . denote string world-sheet coordinates. 

First note that our unperturbed background is under- 
going FRW expansion in the non-compact dimensions 
while the compact direction is stabilized at the self dual 
radius by a string gas. As a preliminary, we wish to show 
that long wavelength perturbations of the metric do not 
affect the motion of strings to first order in the pertur- 
bation variables. We begin with the Polyakov action for 
a closed string 

S = J d 2 ad a X A d a X B g AB (X) , (Al) 

where we work in conformal gauge from the outset. Con- 
sider now perturbations around the background: qar = 
9ab + ^ab- It was shown in the appendix of [42] how 
one can inherit the string spectrum and the constraint al- 
gebra provided the derivatives of the background metric 
9ab (i- e - the Hubble factors) are several orders of mag- 
nitude smaller than the string energy scales. We want to 
make sure that we stay within this regime while studying 
the dynamical compactification of the extra dimension. 
Thus, when we consider the equations of motion for a 
closed string in this perturbed spacetime, we must ne- 
glect all derivatives of the unperturbed metric compared 
to the worldsheet derivatives of the string co-ordinates, 
as these will be of the order of the string scale whereas 
the metric derivatives are constrained to be much smaller 

H. 

Consider the equation of motion for a closed string 

d a (d a x A g CA (x)) = ^d a x A d a x B d c gA B (x), 

from which follows 

x -d a x A d a x B d c h A B = d a d a x c + h c A d a d a x A 

+d a x A d a x B d A h B DgE c ■ 

Now we expand X A as X A = X A + 5X A , where X A 
is the solution to the unperturbed equation of motion 
d a d a X A — 0, see appendix of 01 • Our perturbed equa- 
tion of motion then becomes 

^d a x A d a x B d c h AB = d a d a sx c (A2) 

+d a X A 8 a X B d A h B D9% C ■ 

Notice that there are two vastly different scales in the 
equation above: the worldsheet derivatives and the 
spacetime derivatives. If we imagine expanding Ylab in 
terms of fourier modes, the derivatives would bring down 
a factor of the 4-d wave- vector. Realizing that we need to 
stay within the domain of validity of our setup, namely 
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that we cannot introduce fluctuations in the metric that 
would invalidate the conclusion in the appendix of |42j 
concerning the string spectrum, we need these fluctua- 
tions to vanish on energy scales several orders of magni- 
tude smaller than the string scale. This still allows for 
metric fluctuations that are considerably short distance 
in nature, but nevertheless above the string scale. With 
this in mind, we can neglect all the terms in (| A2|) that 
are multiplied by derivatives of the metric perturbation, 
since these will be orders of magnitude smaller than the 
other ones. Henceforth, we obtain the result 

d a d a 8X A = , (A3) 

that is, the string motion is unchanged by long wave- 
length perturbations. 



Ordinarily this would signal the end of any hopes to 
study perturbative physics; this is not the case for us, 
since we are primarily interested in perturbing around a 
background which consists of a gas of such strings. As 
a result, if we consider the perturbed energy- momentum 
tensor for such a gas, it will contain no anisotropic stress. 
To see this, consider the off-diagonal spatial components 
which can be written as 

ST) « {p l 5 Pl ) + (6p l Pj ) , (A4) 

where i ^ j and (...) indicates an ensemble (gas) average. 
These components vanish, since we just determined that 
for individual strings 8p l = 0, which means that there is 
no anisotropic stress to deal with in our setup. 
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